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INTRODUCTION 


Traditionally, long term measurements of atmospherically propagated sound signals 
have consisted of time series of multiminute averages. Only recently have continuous 
measurements with temporal resolution corresponding to turbulent time scales been 
available. With modern digital data acquisition systems we now have the capability to 
simultaneously record both acoustical and meteorological parameters with sufficient 
temporal resolution to allow us to examine in detail relationships between fluctuating sound 
and the meteorological variables, particularly wind and temperature, which locally 
determine the acoustic refractive index. 


The atmospheric acoustic propagation medium can be treated as a nonlinear 
dynamical system, a kind of signal processor whose innards depend on thermodynamic 
and turbulent processes in the atmosphere. The atmosphere is an inherently nonlinear 
dynamical system. In fact one simple model of atmospheric convection, the Lorenz 
system^), may well be the most widely studied of all dynamical systems. In this paper we 
report some results of our having applied methods used to characterize nonlinear dynamical 
systems to study the characteristics of acoustical signals propagated through the 
atmosphere. For example, we investigate whether or not it is possible to parameterize 
signal fluctuations in terms of fractal dimensions. For time series one such parameter is the 
limit capacity dimension. Nicolis and Nicolis were among the first to use the kind of 
methods we have to study the properties of low dimension global attractors^ 2 ). 

In this paper we show, for example, that the limit capacity dimensions for 
atmospherically propagated acoustic signals are greater than those of either the wind speed 
or the along (propagation) path wind component. Turbulence is the phenomenon which 
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most strongly controls fluctuations in the acoustic refractive index t| . Variations in acoustic 
refractive index are a function of velocity, temperature and, to a lesser extent, humidity 
fluctuations. Written in terms of the turbulent structure function parameters and neglecting 
humidity, variations in r) are 

C* c =aC* + pC^+7C£C> (1) 

where a and p are constants. Gamma is not a constant but rather a function, in particular, 
of the stability (heat flux). 

Although the use of nonlinear dynamical methods is now rapidly growing, they are 
not yet nearly so widely known as, e.g., linear Fourier methods^). Thus we summarize 
here the basic analysis method as well as the results of using it. 


DIMENSIONAL ANALYSIS OF A TIME SERIES: SOME FUNDAMENTALS 


When one is working out of doors it is virtually impossible to measure all of the 
potentially important environmental variables. Nevertheless it may be possible to 
extract most of the information necessary to define signal variability by analyzing 
appropriately combined acoustic and meteorological measurements. 

Takens' theorem ( 4 ) defines the largest embedding dimension which is needed 
to analyze a single time series and, thus, to obtain an accurate fractal dimension for the 
system. The embedding dimension is the state space in which an object can be 
visualized. For any system having a fractal dimension, e.g. the well known Lorenz 
attractor, Takens' theorem states that a maximum embedding dimension of 2d+l is 
needed, where d is the fractal dimension rounded to the next higher integer. Thus an 
embedding dimension of seven should define the Lorenz system, which has a fractal 
dimension of 2.06. A system might be described in fewer dimensions, but Takens' 
theorem sets an upper bound for the state space in which the attractor can be embedded. 

In the analysis of a time series, if an embedding dimension is used which is less 
than prescribed by Takens' theorem, the fractal dimension may not be saturated (i.e., 
reached its peak value). However, as schematically shown in figure 1, if an embedding 
dimension of higher order is used, little, if any additional information will be gained^ 5 ). 


Practically it is important to work with the minimum required embedding dimension in 
order to minimize computational costs. 



Figure 1. Fractal dimension as a function of embedding dimension 


Lagging 


In order to extract all information contained within a single time series it is 
necessary to reconstruct m single order equations. Let the single time series X=F(t) be 
the set of points (xi,x 2 ,x 3 , ...) which are separated by a distance Ax. First we 

approximate the first derivative of F(t) to be 

X ' = = F ' +1 ~ F ‘ ( 2 ) 

dx Ax 

Actually, there is redundant information in the first derivative, as Fj is the original time 

series. Therefore, an embedding dimension of two space is created when the original 
time series is shifted by one time step (Fj, Fj+i)T. For higher order systems this 

process is continued until one has created a state space which is large enough so that the 
attractor can be unfolded. 


If the spacing between points in the approximation of the derivative is too small 
then points will appear to be totally correlated and cannot be considered as independent 
coordinates^). Similarly, if the spacing is too great adjacent points will appear to be 
unrelated (see figure 2). 



In practice, instead of using successive points in the time series to calculate the 
derivatives, the time series is lagged by a certain number of points. Lagging consists of 
setting X' = F'(t) equal to the i+fo sample of F(t), (i.e. i, i+/, i+2 /, ...), where / is 
the size of the lag. Lagging the time series allows one to form a matrix as 


Fi 1 


r x,, x 2 , x 3 , •••"! 


f x l 


F i+ , 
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X |+I> X 2 + l> X 3+J’ *” 
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where each column of the matrix defines a single point in (2d+l) phase space. 



Figure 2. Lorenz attractor with correlated and uncorrelated lags 


To estimate an appropriate value of the lag size, three methods are commonly 
used( 7 ), the autocorrelation time, mutual information, and visualization. For systems 
having an unknown fractal dimension the autocorrelation method appears to be the 
conservative approach. To determine the lag there are two possible ways of 
interpreting a graph of the autocorrelation time scale as shown in figure 3. One is to 
take the point halfway to the first zero crossing. A second approach is to determine the 
halfway point to where the autocorrelation curve becomes parallel to the x-axis. If no 
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zero crossing exists this is the only practical method. For the example shown the two 
methods yield lag sizes of 52 and 35, respectively. 


Another method for estimating a proper lag is called mutual information. In this 
case one increases the lag size until no new information is gained and then defines that 
point as being the appropriate lag size. 


Autocorrelation 



Figure 3. Autocorrelation of the Lorenz attractor 


What is called visualization or visual reconstruction may also be used. 
Visualization is often used in situations where one has some prior knowledge of the 
system. This method consists of graphically reconstructing the attractor with various 
lags. If the topology of the attractor is known, e.g., as in the case of the Lorenz 
attractor, then the lag that appears closest to that for the real system is determined to be 
the appropriate one. 


Limit Capacity Dimension 


The limit capacity is one of four commonly used fractal dimensions: capacity, 
correlation, information, and Lyapunov. Determination of the limit capacity dimension is 



made as follows^). If one lets N(e) represent the minimum number of m-dimensional 
cubes of length 8 needed to enclose the time series, then as 8 decreases one expects N(b) to 
increase. 


N(s) oc s d< 


and therefore the capacity dimension is defined as: 



= lim 

£— >0 


log[AT(g)] 

io 4%. 


(5) 


The output of this limit capacity algorithm gives a lower bound to the dimension of the 
attractor. 


Determining the Dimension from Graph 


Figure 4 shows the result of applying the limit capacity algorithm to the Lorenz 
system. In this representation the bin numbers represent distances between points on the 
attractor. To estimate the limit capacity dimension of a data series a "stable plateau" region 
must be determined. Definition of the stable plateau may be subjective. Definition of 
optimal methods for determining the dimension and hence uncertainties in it are still being 

researched^). 



Figure 4. Limit capacity dimension of Lorenz attractor 
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ACOUSTICAL AND WIND SIGNALS STUDIED 


We have applied the same methods as described above for the Lorenz system to 
hour long recordings of constant frequency sound signals, and to wind speed and the along 
path component of the wind. These measurements were made as part of the comprehensive 
Joint Acoustic Propagation Experiment (JAPE) study. Acoustic receivers were located 1 
Km from the sound source and spaced logarithmically on a tower to a height of 32 meters. 
Three tones of 80, 200, and 500 Hz were transmitted. For this analysis the original 2048 
samples per second were averaged to one quarter second. The corresponding wind time 
series had one-tenth second resolution. So far dimensional calculations for only the 80 and 
500 Hz tones have been completed. 


Lags were determined by calculating the autocorrelation time for each time series. 
Appropriate lags for the acoustic transmission loss (TL) signals varied between 12.5 and 
187.5 seconds; lags for the wind signals ranged from 150 to 400 seconds. Wind speeds 
were less than 6 meters per second during the recording period. Correlation times of the 
TL signals measured at 0, 2, and 32 meters decreased both with height and frequency 
(figure 5). 



Figure 5 . Correlation times vs. height for acoustic signals 




Recall that for the Lorenz attractor an embedding dimension of seven was sufficient 
to unfold the attractor. However, for our acoustic and wind signals the dimension of the 
attractor was unknown. Thus it was necessary for us to calculate the limit capacity for a 
number of different embedding dimensions (figure 6). This was done repeatedly until it 
appeared that the dimensional information had saturated. Figure 6 shows the progression 
of the limit capacity dimension with increasing embedding dimensions until saturation was 
reached at roughly an embedding dimension of 12. 



Figure 6: Dimension of 80 Hz tone at 2 meters for various embedding dimensions 


Results 


Tables I and II, respectively, summarize the calculated limit capacity dimensions for 
the acoustic signals, the wind speed and the along path component of the wind. 


114 





* higher embedding dimensions were not able to be used due to the limited data set and 
high lag. 


Graphical representation of the change in limit capacity dimension with height is 
shown in figure 7 . 
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Limit capacity dimension of acoustic and wind speed signals 



Figure 7. Limit Capacity of acoustic and wind signals 


CONCLUSIONS 


Low order limit capacity dimensions have been determined to exist for both the 
acoustic and wind time series. These results confirm the existence of local attractors. The 
acoustical multivariable dependent signals have higher order attractors than were found for 
the independent meteorological input variables. 


The limit capacity dimension of the acoustic signals appears to increase with height 
and frequency. We believe that this is due to the role which large eddies (thermals) in the 
convective boundary layer (CBL) play in controlling intermittent space-time variations in 
the acoustic refractive index. The properties of propagated sound are sufficiently sensitive 
to those eddies so that tomographic methods may be used to indirectly measure their 
properties (10). 
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We expect that with further dimensional analysis it will be possible to define low 
order dynamical models that will more precisely define the variability of acoustic signal 
fluctuations than can be done presently with linear methods. Further studies will require, 
however, several multihour time series recorded in both stable and unstable boundary layer 
conditions. The single hour time series recorded during JAPE is of insufficient length. 
Since large eddies appear to be the dominant signal controlling mechanism, it would also be 
helpful to have measurements over transmission paths ranging from about 2.5 to 10 km. 
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